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A supermassive black hole can disrupt a star when its tidal field exceeds the star's self-gravity, 
and can directly capture stars that cross its event horizon. For black holes with mass M > 1O 7 M , 
tidal disruption of main-sequence stars occurs close enough to the event horizon that a Newtonian 
treatment of the tidal field is no longer valid. The fraction of stars that are directly captured is 
also no longer negligible. We calculate generically oriented stellar orbits in the Kerr metric, and 
evaluate the relativistic tidal tensor at the pericenter for those stars not directly captured by the 
black hole. We combine this relativistic analysis with previous calculations of how these orbits are 
populated to determine tidal-disruption rates for spinning black holes. We find, consistent with 
previous results, that black-hole spin increases the upper limit on the mass of a black hole capable 
of tidally disrupting solarlike stars to ~ 7 x 10 8 Mq. More quantitatively, we find that direct stellar 
capture reduces tidal-disruption rates by a factor ~ 2/3 (1/10) at M ~ 10 7 (10 8 )Mq. The strong 
dependence of tidal-disruption rates on black-hole spin for M > 10 8 Mq implies that future surveys 
like the Large Synoptic Survey Telescope that discover thousands of tidal-disruption events can 
constrain supermassive black-hole spin demographics. 



I. INTRODUCTION 

In 1943, active galactic nuclei (AGN) were discovered 
with emission lines Doppler broadened to widths > 1, 000 
km/s pp. Twenty years later, theorists proposed that 
these AGN were powered by accretion onto compact ob- 
jects of masses 10 5 — 1O 8 M [2]. Such massive objects 
cannot support themselves against gravitational collapse 
into supermassive black holes (SBHs) [3]. SBH masses 
are tightly correlated with the luminosity 4J, mass [5], 
and velocity dispersion [BJ of the spheroidal component 
of their host galaxies. 

SBHs primarily grow by accreting gas driven into 
galactic centers by tidal torques during major mergers 
[3 [8]. However, SBHs can also grow by directly captur- 
ing stars that cross their event horizons or by accreting 
debris from stars passing close enough to be tidally dis- 
rupted 0. Such tidal-disruption events (TDEs) could 
also power bright flares of radiation as the stellar de- 
bris is shock heated and subsequently accreted [THl E] ■ 
Several potential TDEs have been found in x-ray [T2] . 
UV [T3J, and optical [TJj surveys, and tidal debris may 
also fuel recent blazar activity seen by the Swift satellite 
[ToT - fTS] . The handful of TDEs found in current surveys 
implies that thousands more may be found each year in 
future surveys by the Large Synoptic Survey Telescope 
fL4l 119] . A detailed comparison between predicted and 
observed TDE rates will provide important constraints 
on SBHs and the galactic centers in which they reside. 

Frank and Rees [S] were among the first to estimate 
TDE rates. They introduced the concept of a "loss cone" 
in the stellar phase space that would be depopulated 
by tidal disruption within a dynamical time tdyn- Stars 
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within the loss cone have velocities lying in a cone about 
the radial direction. TDE rates are set by the rate at 
which stellar diffusion from other portions of phase space 
refills this loss cone. Frank and Rees evaluated stellar 
fluxes into the loss cone at a critical radius r cr i t at which 
stellar diffusion operating on a reference time scale tn [20] 
could refill the loss cone within idyn- Cohn and Kulsrud 
[2"T] provided a more sophisticated treatment of stellar 
diffusion into the loss cone by numerically integrating 
the Fokker-Planck equation in energy- angular momen- 
tum space. More recently, Wang and Merritt [55] have 
revised predicted TDE rates using more realistic galactic 
density profiles and the observed relation between SBH 
mass and host-galaxy velocity dispersion [5J. 

These analyses focused on smaller SBHs for which a 
Newtonian treatment of tidal forces is valid and the num- 
ber of directly captured stars is negligible compared to 
the number that are tidally disrupted. Manasse and Mis- 
ner [23J introduced Fermi normal coordinates that are 
ideal for a relativistic treatment of the tidal tensor, and 
calculated this tidal tensor for radial geodesies of the 
Schwarzschild metric for nonspinning SBHs. Marck |24j 
generalized this calculation to generically oriented time- 
like geodesies of the Kerr metric [25] for spinning SBHs. 
Beloborodov et at [5B] used this tidal tensor to calcu- 
late the relativistic cross sections for tidal disruption for 
a range of initial orbital inclinations with respect to the 
SBH spin. Ivanov and Chernyakova 27J used a numeri- 
cally fast Lagrangian model of a tidally disrupted star to 
investigate how stellar hydrodynamics affects these rela- 
tivistic cross sections. In this paper, we combine a sim- 
ilar relativistic treatment of tidal disruption and direct 
capture with existing calculations of loss-cone physics to 
derive improved predictions of TDE rates for massive and 
highly spinning SBHs. 

The first step in calculating TDE rates is to estab- 
lish criteria for determining when a star is tidally dis- 
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rupted. Most tidally disrupted stars are initially on 
highly eccentric or unbound orbits characterized by the 
distance r of their pericenters from the SBH. An order- 
of-magnitude estimate of the maximum value of r for 
which tidal disruption occurs can be obtained by equat- 
ing the differential acceleration GMR^/r 3, experienced 
by a star of mass and radius i?* in the tidal field 
of a black hole of mass M to the star's self-gravity 
Gm*/R%. This implies that a star will be tidally dis- 
rupted when r < ttd — (M/m*) 1 / 3 ^*. A star will be 
directly captured by the SBH when r is less than the ra- 
dius of the event horizon, which for a nonspinning SBH 
is equal to the Schwarzschild radius rg = 2GM/c 2 . Since 
ttd oc A/ 1 / 3 while rs oc M, the ratio of tidally disrupted 
to directly captured stars will decrease with increasing 
SBH mass. Equating these two distances, we find that a 
SBH with mass M greater than 
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should directly capture stars instead of tidally disrupting 
them. 

Our estimate of ttd assumed that the gravitational 
field of the SBH was that of a Newtonian point par- 
ticle, which should only be valid for t^d 3> r s- One 
should be very suspicious of using this estimate at the 
event horizon, as we did when deriving M max above. In 
a proper general-relativistic treatment, the spacetime of 
a spinning SBH is described by the Kerr metric |25j . 
The Kerr metric is a two-parameter family of solutions 
to Einstein's equation fully specified by the SBH mass 
M and dimensionless spin a/M < 1. Theoretical esti- 
mates of SBH spins depend sensitively on the extent to 
which SBHs grow by accretion or mergers. SBHs accret- 
ing from a standard thin disk |28j can attain a limiting 
spin a/M ~ 0.998 29 after increasing their mass by a 
factor - a/6 [30]. The spins of SBHs formed in mergers 
vary greatly depending on whether the spins of the ini- 
tial binary black holes become aligned with their orbital 
angular momentum prior to merger |31| . SBH spins can 
be inferred from observations of iron Ken lines in AGN 
x-ray spectra |32j . Large spins have been measured, such 
as a/M = 0.989;to;oo2 in the Seyfcrt 1.2 galaxy MCG-06- 
30-15 [33], although reliable estimates are only available 
for a small number of systems. 

Given the large sample of observed TDEs expected in 
the near future and the wide range of predicted SBH 
spins, it is important to determine the extent to which 
TDE rates depend on SBH spin. This is the primary goal 
of this paper. The greater the spin dependence, the more 
tightly observed TDEs will be able to constrain the distri- 
bution of SBH spins. In Sec. [Hj we review how the tidal 
field is calculated along timelike geodesies of the Kerr 
metric. In Sec. |III[ we describe the Monte Carlo simula- 
tions we performed to determine which stellar orbits lead 



to tidal disruption. We then use these simulations to cal- 
culate expected TDE rates in Sec. |IV| The implications 
of our findings are discussed in Sec. |V| 



II. TIDAL FIELDS ALONG KERR GEODESICS 

In Boyer-Lindquist coordinates |34j and units where 
G = c = 1, the Kerr metric takes the form 
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where £ = r 2 + a 2 cos 2 9 and A = r 2 — 2M r + a 2 . This 
metric is both stationary (independent of t) and axisym- 
metric (independent of 0). Massive test particles travel 
on timelike geodesies of the Kerr metric. Individual stars 
have masses ra* ~ M much less than those of SBHs 
(1O 6 M < M < 1O 1O M ), and radii i?* - R Q ~ 7 x 10 10 
cm less than the Schwarzschild radius 
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We can therefore consider them to be test particles for 
the purpose of determining their orbits. The position 
(r, 9, </>) of a star as a function of proper time t evolves 
according to the equations 35 
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where the specific energy E, angular momentum L Z1 and 
Carter constant Q are conserved along geodesies. 

Although Boyer-Lindquist coordinates reduce to flat- 
space spherical coordinates in the limit r —¥ oo, the 
nonzero off-diagonal elements of the Kerr metric ([2]) im- 
ply that these coordinate vectors do not constitute an 
orthogonal tetrad at finite r. The gravitational-field gra- 
dients experienced by freely falling observers are more 
conveniently expressed by projecting them onto an or- 
thonormal tetrad A M like that provided by Fermi normal 
coordinates [23]- This coordinate system can be used 
to specify points in the neighborhood of a central time- 
like geodesic, such as that traversed by a star orbiting 
a Kerr SBH. The timelike member of this tetrad Ao is 
the tangent vector along the central geodesic, while the 
spacelike vectors A^ (i = 1, 2, 3) span the plane in the tan- 
gent space othogonal to Ao- The point (r, x l ) in Fermi 
normal coordinates is reached by starting at the location 
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of the star at proper time r and moving a proper dis- 
tance R = \/J2i( xl ) 2 along the spacelike geodesic whose 
tangent vector is J2i x% \- 

In Fermi normal coordinates, the time-time component 
of the metric can be Taylor expanded as 



where 



(5) 



where R a p-yS is the Riemann curvature tensor projected 
onto the orthonormal tetrad A M . This implies that the 
tidal potential ^tidai experienced by a star is 



$tidai - -\{g^K + 1) = \ r > j ■'■'■>■' + ■•■ . (6) 



where Cy = Rmoj is the tidal tensor. Although the 
higher-order corrections to the tidal potential denoted 
by the ellipsis can sometimes be significant [36] . in this 
paper we consider only the term quadratic in x l . The 
tidal tensor Cy is a symmetric, traceless 3x3 matrix 
whose eigenvectors denote the principal axes along which 
the star is stretched or squeezed, and whose eigenvalues 
denote the extent of this stretching and squeezing. 

The problem of calculating the tidal potential $tidai 
thus reduces to choosing an appropriate orthonormal 
tetrad A M for generic Kerr geodesies and projecting the 
Riemann tensor onto this tetrad. This has already been 
accomplished for us by Marck [53], who found 
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The angle ^ evolves along the geodesic to ensure that Ai 
and A 3 are parallel propagated. 

The fully general tidal tensor of Eq. Q is intimidating, 
but we can gain insight by considering the tidal tensor 
for equatorial geodesies (9 — tt/2, Q = 0) whose nonzero 
elements are 



1 - 3- 



C u = 

Cl3 = —3 -g 



K 



cos 2 * 



f 
2 + K 



M 



-M cos* sin* 



C22 = 

C33 = 

The eigenvalues 
3K/r 2 )M/r 3 , and 
force is 

F = 



K\ M 



r 2 +K 



sin 2 * 



M 



(9a) 

(9b) 
(9c) 

(9d) 



of this tensor are 
-2(l + 3K/2r 2 )M/r 3 . 



tidal 



Cij 



M/r 3 , (1 + 
Since the tidal 



(10) 



the positive eigenvalues correspond to directions in which 
the star is squeezed while the negative eigenvalues cor- 
respond to the direction in which it is streched. In the 
Newtonian limit K/r 2 — > 0, the eigenvalues reduce to 
—2M/r 3 and the doubly degenerate eigenvalue M/r 3 . 
This degeneracy reflects the restoration of symmetry be- 
tween the 9 and <f> directions at large r, where the effects 
of the SBH's spin are negligible. Stretching occurs in the 
radial direction corresponding to the eigenvalue —2M/r 3 . 
Note that despite one's possible intuition to the contrary, 
the tidal force remains finite at both the innermost stable 
circular orbit and even the event horizon itself. 

To determine whether a star on a given orbit is 
tidally disrupted, we check at the pericenter of that orbit 
whether the outward tidal force in the direction corre- 
sponding to the negative eigenvalue of the tidal tensor 
exceeds the inwards Newtonian self-gravity of the star. 
We assume that the tidal field is maximized at the peri- 
center as in the Newtonian limit. If j3— denotes the nu- 
merical value of this eigenvalue, tidal disruption occurs 
if 
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In the Newtonian limit /?_ = — 2M/r 3 discussed above, 
this condition is equivalent to the more familiar expres- 
sion 
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Although our condition ( 11 ) for tidal disruption is only 



approximate, we expect it to be conservative for several 
reasons. It neglects that the tidal force has already raised 
bulges on the star's surface before the star reaches the 
pericenter, so the stellar radius i?* appearing in Eq. ( 11 ) 



should exceed its value in hydrostatic equilibrium far 
from the SBH. It also assumes that the star is nonro- 
tating, while in reality the torques exerted on the tidally 
distorted star will cause it to partially corotate with its 
orbit. These torques are likely to be complicated for 
a generic nonequatorial Kerr geodesic, but we can gain 
some insight by again considering the Newtonian limit. 
Stars rotating with angular velocity Q will be disrupted 
at a radius 
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in this limit. For corotating stars on circular orbits (f2 2 = 
GM/r 3 ), the first factor in parentheses on the right-hand 
side of Eq. ( 13 ) equals 3 as in the definition of the radius 
of the Hill's sphere [37]. For a star corotating at the 
pericenter of a parabolic orbit like that expected for a 
star approaching a SBH from a large distance, this factor 
equals 4. Our assumption that the star is nonrotating is 
conservative because the condition ( 13 ) is most restrictive 
for — 0, although ttd(^) only varies by the modest 
factor 2 1 / 3 . 



Our criterion ( 11 ) might overestimate the rate at which 



stars are fully disrupted since they might lose their outer 
layers while maintaining their dense cores. In the New- 
tonian limit, Phinney |38| showed that stars will not be 
fully disrupted until 
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where k is the constant of apsidal motion and fGm* /R* 
is the star's binding energy. The factor k/f = 0.3(0.02) 
for stars with convective (radiative) atmospheres, but the 
exponent of 1/6 ensures that ttd is only weakly depen- 
dent on this factor. We ignore this factor and keep our 
criterion (11) for the remainder of this paper, but the 



Monte Carlo simulations described in the next section 
could easily be evaluated with a new criterion that incor- 
porates this factor or a different choice of stellar proper- 
ties than to* = Mq, i?* = R & . 

A detailed study of the fraction of mass loss as a func- 
tion of SBH and orbital parameters is beyond the scope 
of this paper, but such a study for selected orbital incli- 
nations has been performed by Ivanov and Chernyakova 
|27j . They recognized that for a given SBH mass and 
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FIG. 1: The mass Af max of the heaviest SBH capable of dis- 
rupting a star of solar mass and radius as a function of SBH 
spin a/M. The red squares show the values listed in Table 
2 of |27| derived using a simple hydrodynamical model. The 
solid blue curve sho ws our prediction according to the rela- 
tivistic criterion (111, while the dashed black curve shows the 



Newtonian prediction ( 12 1 



spin, stars on prograde, equatorial, marginally bound or- 
bits are the most likely to be tidally disrupted. Using 
a simple but computationally inexpensive hydrodynam- 
ical model, they calculated the mass M max of the heav- 
iest SBH capable of tidally disrupting stars without di- 
rectlycapturing them. As in the Newtonian prediction of 
Eq. |l|, M max oc to* 1 ^ 2 R'J 2 . In Fig.Jl] we compare their 
predictions (red squares) to our own using the relativistic 
criterion ( |11[ ) (solid blue curve) and Newtonian criterion 
( 12 ) (dashed black curve) for stellar mass to* = Mq and 
radius _R* = Rq. We see that the relativistic correc- 
tion to the Newtonian prediction is significant, and that 
our simple criterion ( 11 ) does a reasonable job given the 



~ 50% uncertainty in the simulations [2"T| . 

We see that in the maximally spinning limit (a/M — > 
1), a SBH as massive as ~ 1O 9 M is capable of tidally 
disrupting main-sequence stars. This prediction is consis- 
tent with earlier simulations [301 SI] that demonstrated 
this possibility. The above scaling of M max with stel- 
lar mass and radius suggests that a white dwarf with 
to* ~ Mq,R* ~ O.Oli? could be tidally disrupted by 
a maximally spinning SBH as massive as 10 6 M Q . This 
conclusion helps alleviate tension between the small SBH 
mass required for the interpretation of Swift J1644+57 as 
a white-dwarf tidal disruption [42] and the larger value 
of M inferred from the relation between SBH mass and 
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the equations of motion Q with the star located at an 
initial position (r, 9, <p) in Boyer-Lindquist coordinates. 
Since these equations are independent of 4> as is the tidal 



FIG. 2: Our choice of coordinates for determining the ini- 
tial conditions for integrating the equations of motion Q for 
stellar orbits. The SBH is located at the origin, and the star 
is located at Boyer-Lindquist coordinates (r,9,(f>). <d is the 
angle between the stellar velocity v and the inwards radial 
direction — f , while $ is the angle between the component of 
v perpendicular to f and the unit vector eg in the 6 direction. 



host-galaxy velocity dispersion [6 . 



III. MONTE CARLO SIMULATIONS 

Unlike the Newtonian two-body problem, there is no 
general analytic solution to the relativistic equations of 
motion Q. We must integrate these equations explicitly 
for every orbit we consider. Stars that will eventually be 
tidally disrupted are scattered onto their final orbits at 
radii r ^> ttd- These orbits may or may not be grav- 
itationally bound to the SBH, but their Newtonian or- 
bital energies ~ m^a 2 , where a is a typical velocity at 
r 3> rxD, are much less than the rest-mass energy m^c 2 . 
It is therefore an excellent approximation to set the spe- 
cific energy E appearing in Eqs. (El equal to unity in 
units where c = 1. The Kerr metric (21) is axisymmetric, 
so our results are independent of the initial value of <f>. We 
must perform Monte Carlo simulations with an appropri- 
ate distribution of the remaining variables {r,9, L Z ,Q} 
to determine what fraction of orbits are tidally disrupted 
according to our relativistic criterion (11), where the neg- 



ative eigenvalue /?_ of the tidal tensor Cy depends on all 
these variables. 

We illustrate the geometry of the problem and our 
choice of coordinates in Fig. [2j We begin integrating 



tensor Cy , we do not actually need to integrate Eq. ( 4c 



We choose an initial radius r = 2000M, where relativisic 
corrections are small, and check that our results are in- 
sensitive to this choice. In this limit, the constants of 
motion are given by 
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where v — (2M/r) 1 ' 2 is the magnitude of the initial stel- 
lar velocity and 9 and $ are the angles described in 
Fig. [2j Since the stars at r > ttd do not know about 
the direction of the SBH spin, the stellar distribution 
is axisymmetric about f and there is a uniform distri- 
bution in $. Although astrophysical spheroids do not 
necessarily have isotropic velocity dispersions at large r, 
the stars approaching this close to the SBH belong to the 
tiny fraction of the distribution where the velocity lies in 
a loss cone centered about the radial direction [5]. Since 
there is no reason to expect the distribution function to 
be varying strongly in this small portion of phase space, 
there is a uniform distribution in — 1 < cos 9 < 1 as well. 
However, since the rate at which stars of velocity v enter 
the sphere of radius r is proportional to v • f , we weight 
our distribution by cos 9 during our Monte Carlo sim- 
ulations of stellar orbits. We choose a maximum value 
9 max to avoid wasting computational time on orbits that 
do not closely approach the SBH. 

With this choice of initial conditions, we integrate the 
equations of motion Q until the star reaches the percen- 
ter. We then calculate and tabulate the negative eigen- 
value /3_ of the tidal tensor Cy (JtJ) . We also tabulate 
which stars are directly captured by the SBH when their 
orbits encounter the SBH's event horizon. We integrated 
250,000 stellar orbits for each of several SBH spins, with 
an additional 250,000 with a smaller choice of 9 max to 
increase our sampling of the small number of orbits that 
lead to tidal disruption when M — > M max . 



IV. TDE RATES 

Given a stellar phase-space distribution function, it is 
reasonably straightforward to calculate the rate of TDEs 
using the Monte Carlo simulations described in the pre- 
vious section. In Sec. |IV A| below, we calculate the TDE 
rate as a function of SBH mass assuming that the stars 
approach a Maxwellian distribution with fixed number 
density and velocity dispersion far from the SBH. This 
calculation illustrates the dependence of TDE rates on 
SBH spin. However, astrophysical SBHs reside in galaxy 
spheroids whose properties are tightly correlated with 
SBH mass [2HS]- In addition, the very luminous early- 
type galaxies that host the most massive SBHs have cored 
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profiles at their centers unlike the power-law profiles that 
characterize less luminous early-type galaxies and late- 
type bulges 03]. Predicted TDE rates are sensitive to 
whether galactic centers arc described by cored or power- 
law profiles [25]. Recent observations [H] suggest that 
even the nuclear star cluster at our own Galactic center, 
long believed to have a cuspy profile (p oc r -7 / 4 ) |3S], 
may in fact have a core of radius r CO r G — 0.5 pc |46) . 
Given these uncertainties, it is difficult to make precise 
estimates of astrophysical TDE rates. Despite this, in 
Sec. |IVB| we calculate the TDE rate assuming galaxies 
have isothermal (p oc r~ 2 ) profiles at their centers and 
host SBHs with masses correlated with their velocity dis- 
persions. The results of this calculation shown in Fig. [4] 
illustrate how SBH spin affects TDE rates. 



the "Newtonian" prediction would be that a star is di- 
rectly captured by the SBH if the pericenter of its or- 
bit is less than the Schwarzschild radius pi). The peri- 
center of a parabolic (En — 0) orbit with specific an- 
gular momentum Ln is L N /2GM. Equating this to 
the Schwarzschild radius, a star is directly captured if 
L N 

plies a capture rate 



< L cap = 2GM/c, which according to Eq. (201 im- 



r = 

1 cap 



dT _ (32n)^ 2 n(GM) 2 

— — dL N — 



= 2.1 x 10~ 



M 



yr 



(J 

100 km/s 



V 10 6 M Q 

-l 



10 5 pc~ 



(21) 



A. Maxwellian distribution 

Assume that stars far from the SBH have a Maxwellian 
distribution function 
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with number density n and velocity dispersion a. The 
differential rate at which stars with Newtonian specific 
energy En and angular momentum L n enter a sphere of 
radius r is given by 
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where the volume integral extends over the region v z > 
and E' and L' are given by 
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A star will be tidally disrupted if the pericenter of its 
orbit is less than the tidal-disruption radius r^D (12 1, 
which implies an angular momentum 
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This implies a TDE rate 
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This rate agrees with that in Eq. (16b) of Frank and Rees 
[S] which applies when the critical radius at which the 
loss cone refills on a dynamical time exceeds the SBH's 
radius of influence. If TDEs can only be observed when 
the tidal debris is not directly captured by the SBH, the 



We can use Eq. ( 16 ) and the delta functions to evaluate observed TDE rate will be T T d = TJ, 



the integral to find 
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cap- Since 

r TD oc M 4 / 3 while r cap oc A/ 2 , the TDE rate will vanish 
for M > M max at which ttd = 

We can use this same differential rate dT / 8Ln ( 20 1 to 



If cr <C c, orbits near the SBH will be insensitive to the 
value of En and we can integrate over this variable to 
yield a differential rate 
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calculate the relativistic direct-capture and TDE rates. 
However, we must now rely on the Monte Carlo simula- 
tions of Sec. |III| to determine whether stars are directly 
captured or tidally disrupted, instead of the simple New- 
tonian expressions for L cap and Ltd given above. The 
simulated orbits have a maximum angular momentum 
i max = (2GMr) 4 / 2 sin O max . The total rate at which 
stars on these orbits enter a sphere of radius r = 2000Af 



The divergence of this rate as a — > results from gravi- 
tational focusing, which would channel all stars into the 
SBH in the absence of tangential velocities. 

Before proceeding to the relativistic calculation, let 
us review the Newtonian predictions. Although the 
event horizon is fundamentally a relativistic concept, 



Ttnt — 



Jo 



dT 



dLN = 



(8n)^ 2 nGMr . 2 ^ 
- — sur e r 



( 24 ) 

The rate r cap at which stars are directly captured by 
the SBH is found by multiplying this total rate r to t by 
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FIG. 3: The rates at which stars are directly captured (dotted 
lines) and tidally disrupted (solid and dashed curves) by SBHs 
of mass M in constant-density cores with n = 10 5 pc~ 3 and 
<r = 100 km/s. The black curves show the Newtonian rates of 
Eqs. ( 21 1 and ( 23 1, while the colored curves show the relativis- 
tic rates for SBHs with spins a/M = (red), 0.5 (orange), 0.9 
(green), 0.99 (blue), and 0.999 (purple). The capture rates 
mildly decrease with SBH spin, while for M > 10 8 Mq the 
TDE rates greatly increase with SBH spin. 



the fraction F cap of simulated geodesies that cross the 
event horizon. The TDE rate Ttd is similarly found by 
multiplying r to t by the fraction Ftd of orbits that vi- 



olate the relativistic criterion (11) for tidal disruption 



If r and max are chosen large enough, these fractions 
F oc (r sin 2 max ) _1 so that the physical rates are inde- 
pendent of our choice of initial conditions. 

In Fig. [3j we show the direct-capture rate r cap and 
TDE rate Ttd as functions of SBH mass M for our fidu- 
cial choices n = 10 5 pc -3 and a = 100 km/s. The New- 
tonian prediction for r cap underestimates the true rela- 
tivistic capture rate by about a factor of 4. If we had 
used the true specific angular momentum L z — AGM/c 
for marginally bound geodesies of a Schwarzschild SBH 
as the upper limit of the integral in our Newtonian pre- 
diction ( 21 ), we could have nearly reproduced the correct 



relativistic result. The capture rate is nearly independent 
of SBH spin as indicated by the colored dotted lines ly- 
ing almost on top of each other. This is surprising, since 
the specific angular momentum L z for prograde (retro- 
grade) marginally bound equatorial orbits varies from 
4M (-4M) to 2M (-4.828M) as a/M increases from 
to 1. Near perfect cancellation over orbital orientation 
(6, 0, $) must occur for the capture rate r cap to be so 



mildly dependent on spin, but we do not see any obvi- 
ous reason for this to be the case. Young et al. [48] 
calculated the ratio of the capture rate r cap (a) for Kerr 
SBHs of spin a to the capture rate r cap (0) for nonspin- 
ning SBHs. Equation (B2) of their paper shows that this 
ratio is approximately given by 



cap 



(«) 



cap 



(0) 



= 1-0.0820 



-0.0717 



4 

-0.0864 



M 
(25) 

The small numerical values of the coefficients in this 
expression indicate the weak dependence of the direct- 
capture rate on SBH spin; the ratio is between 0.9 and 
unity over the full range of spins < a/M < 1. 

The TDE rate Ttd exhibits a much stronger depen- 
dence on SBH spin, as illustrated by the strongly vary- 
ing solid colored curves in Fig. [3] At small SBH masses, 
ttd ^ M and the TDE rate for all spins converges to the 
Newtonian result as expected. However, as M increases, 
tidal disruption occurs closer to the SBH where the New- 
tonian approximation is increasingly invalid. This is most 
glaringly apparent for masses M > M max of Eq. for 
which tidal disruption would not be possible in the New- 
tonian limit. The true maximum mass, where the solid 
colored curves in Fig. [3] intersect r TD = 0, is given as 
a function of spin in Fig. [I] Since T^d ^ 

M 4 / 3 , these 

massive SBHs are capable of tidally disrupting even more 
stars than their less massive counterparts. Although the 
spins a/M — 0.99 and 0.999 depicted by the blue and 
purple curves in Fig.[3]may seem extreme, the simple sce- 
nario of growing a SBH from a standard thin accretion 
disk leads to a limiting spin a/m ~ 0.998 quite close to 
the purple curve [25]. Although uncertain, cosmological 
predictions for SBH spin distributions can also be peaked 
near these large values [31] . The primary conclusions to 
draw from our analysis are that relativistic corrections to 
the TDE rate can alter predictions by a factor of several 
for M > 1O 7 M , and that they can allow TDEs to occur 
for SBHs as large as ~ 1O 9 M0. 



B. Real galaxies 

Following Frank and Rees [5J, the rates we calculated 
in the preceding subsection assumed that galaxies had 
constant-density cores outside the SBH's radius of influ- 
ence, 



Th = 



™ = 0.43 pc 

G 



M 



(i0 f; i/. J \ 100 km/s 



(26) 



Real galaxies with either power-law or core profiles have 
mass-density profiles p(r) that monotonically decrease 
with radius. This raises the question of what is the ap- 
propriate number density n to insert in our expressions 
for direct-capture and TDE rates. Frank and Rees [9 
argued that the appropriate density to use is that at the 
critical radius r cr j t at which stellar diffusion can refill the 
loss cone of tidally disrupted orbits on a dynamical time. 



A very crude estimate of this density can be made by 
assuming that r cr i t — r^, an assumption roughly true for 
real galaxies as indicated by Fig. 6 of Wang and Mer- 
ritt [52] . If we further assume that the density profile of 
galactic centers is that of a single isothermal sphere, 



2nGr 2 ' 



(27) 



then inserting n = p(rh)/m* into Eq. (23 1 implies 



V 



TD 



1.3 x 10" 



M 



(ro <; \/. 



-2/3 



100 km/s 



(28) 

Wang and Merritt [25] use the isotropic distribution func- 
tion appropriate for a single isothermal sphere to calcu- 
late the true rate at which the loss cone is refilled by 
stellar diffusion. They find that their results are well 
approximated by the fit 



r' TD ~ 2.5 x 10 _3 yr _1 



( 10' 



100 km/s 



7/2 



(29) 

If we combine this estimate with a recent determination 
of the relation between SBH mass and host-galaxy veloc- 
ity dispersion [47] . 



M 



= 7.58 



1O 6 M """V 100km / S 
we arrive at a final TDE rate of 



4.32 



1% ~ 4.8 x lQ- 4 yr _1 



M 



(l() f . 



(30) 



(31) 



in the Newtonian limit. This estimate should be reason- 
able for the power-law galaxies that dominate the total 
TDE rate; the core galaxies that host the most massive 
SBHs have TDE rates r^ D ^ 10 -5 yr -1 about an order 
of magnitude below that of comparable-mass power-law 
galaxies [22] . 

In Fig. [4j we show how the direct capture of stars by 
spinning SBHs changes this prediction. This figure was 
prepared with the same set of Monte Carlo simulations 



described in Sec. Ill Although there are considerable dif- 
ferences between the Newtonian predictions of Eqs. ( 23 ) 
and (31), these differences result from different treat- 



ments of the stellar populations far from the SBH. Wc 
may therefore simply renormalize our relativistic predic- 
tions Ttd — ^TDTtot of the previous subsection by di- 



viding by Eq. (23) and multiplying by Eq. (31) at each 
SBH mass M. Direct capture reduces the predicted TDE 
rate by a factor - 2/3 (1/10) at M = 10 7 (10 8 )M Q . Al- 
though TDEs are very rare for large SBH masses, they 
are still possible for M < M max ~ 10 9 M o . Since SBHs 
with masses M ~ 10 9 M Q predominantly live in galaxies 
with cored profiles, Fig. [4] may somewhat underestimate 
TDE rates at these masses since the stellar density will 
not fall as steeply with r as the single isothermal profile 
of Eq. ([27]. 
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FIG. 4: The rates at which stars are tidally disrupted by 
SBHs of mass M in power-law galaxies obeying the M — a 
relation. The dashed black line is the prediction of Wang and 
Merritt [22 for I^d with an updated M - a relation. The 
colored curves show our relativistic corrections Ftd to this 
prediction. The TDE rate increases with SBH spin, with the 
given curves corresponding to a/M = (red), 0.5 (orange), 
0.9 (green), 0.99 (blue), and 0.999 (purple). 



V. DISCUSSION 

Astronomers have sought to observe the electromag- 
netic flares associated with TDEs ever since this possi- 
bility was proposed by Rees [TU] . Several potential TDEs 
were discovered over the past 15 years by the Roentgen 
Satellite (ROSAT) [12 and the Galaxy Evolution Ex- 
plorer [13], and the recent discovery of additional TDEs 
by both the Sloan Digital Sky Survey [H] and Swift [PoT - 
IT8] has renewed interest in this phenomenon. While in- 
dividual TDEs may provide new insights into SBH accre- 
tion physics, the large samples that may soon be avail- 
able |14| will uniquely probe the whole population of both 
active and quiescent SBHs. While overall TDE rates de- 
pend on stellar populations at galactic centers, the upper 
bound on the mass M of SBHs capable of tidal disruption 
is a sensitive measure of SBH spins. For M > 1O 7 M , 
tidal disruption occurs deep enough in the SBHs poten- 
tial well that Newtonian gravity is no longer valid. Fur- 
thermore, there is no reason to expect the orbital angu- 
lar momenta of tidally disrupted stars to align with SBH 
spins. For both these reasons, accurate calculations of 
TDE rates require evaluation of the relativistic tidal ten- 
sor Cij on a representative sample of generically oriented 
Kerr geodesies. 
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We have performed a series of Monte Carlo simula- 
tions that provide this required sample. We use this 
sample to calculate TDE rates for spinning SBH as a 
function of their mass M, both in constant-density cores 
and in isothermal spheres that approximate real power- 
law galaxies. We find that for M > lO 7 Af , a significant 
fraction of stars will be directly captured by the SBH's 
event horizon instead of being tidally disrupted and sub- 
sequently accreted. This will reduce the observed TDE 
rate assuming that directly captured stellar debris will 
not have the chance to radiate appreciably before being 
swallowed by the SBH. Above M ~ 1O 8 M , only highly 



spinning (a/M > 0.9) SBHs will be able to produce ob- 
servable TDEs. Theory (31] and observation [32l [33] 
suggest that most SBHs may have such large spins, but 
further observations are needed to investigate this pos- 
sibility. A future survey like the Large Synoptic Survey 
Telescope [19j that finds thousands of TDEs may provide 
important constraints on the distribution of SBH spins. 
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